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Abstract: In this letter we announce rigorous results that elucidate the relation between metastable 
states and low-lying eigenvalues in Markov chains in a much more general setting and with considerable 
greater precision as was so far available. This includes a sharp uncertainty principle relating all 
low-lying eigenvalues to mean times of metastable transitions, a relation between the support of 
eigenf unctions and the attractor of a metastable state, and sharp estimates on the convergence of 
probability distribution of the metastable transition times to the exponential distribution. 
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I. Introduction. 

The phenomenon of metastability has been a fascinating topic of non-equilibrium statistical me- 
chanics for a long time. Currently, it has found renewed interest in the investigation of glassy systems 
and aging phenomena which appear to play a central role in many physical and non-physical systems. 
An approach to link metastability to spectral characteristics, in particular low-lying eigenvalues and 
the corresponding eigenfunctions has been proposed by Gaveau and Schulman [GS] . Such an approach 
is appealing not only because it allows to characterize metastability in terms that are intrinsically dy- 
namic and make no a priori reference to geometric concepts such as "free energy landscapes", but 
also since it allows numerical investigations of metastable states via numerical spectral analysis (see 
Schiitte et al. [S,SFHD] for applications to conformational dynamics of biomolecules). 

Relating metastability to spectral characteristics of the Markov generator or transition matrix is 
in fact a rather old topic. First mathematical results go back at least as far as Wentzell [W] (see 
also [M] for more recent results) and Freidlin and Wentzell (see [FW]). Freidlin and Wentzell relate 
the eigenvalues of the transition matrix for Markov processes with exponentially small transition 
probabilities to exit times from "cycles" ; Wentzell has a similar result for the spectral gap in the 
case of certain diffusion processes. All these relations are on the level of logarithmic equivalence, 
i.e. of the form limejo elii(AfT/) = where e is the small parameter, and A-,T/ are the eigenvalues, 
resp. exit times. This rather crude level of precision persists also in the more recent literature and 
prevents, in particular, applications to systems with unbounded numbers of metastable states which 
are particularly relevant for glassy systems. 

In this letter we announce results that - for a large class of Markov chains - improve this situation 
considerably: in particular we allow for the number of metastable states to grow (with e.g. the 'vol- 
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ume'), and we give precise control of error terms for 'finite volume' systems. Moreover, we provide 
representations for all quantities concerned that are computable in terms of certain 'escape probabili- 
ties' that are in turn well controllable through variational representations [BEGKl]. A more detailed 
exposition of our results, as well as the proofs, will be given in two forthcoming papers [BEGK2,EK]. 

2. Stable set and metEistable states. 

We will consider in the sequel Markov chains Xf with state space Fjv, discrete time^ f e N, and 

transition matrices Pjy- We will assiimc that for any fixed they are ergodic, and have a unique 
invariant distribution Qn- We are interested in the situation when these chains exhibit "metastable" 
behaviour; loosely speaking, this means that the state space Gn can be decomposed into subsets SN,i 
such that the typical times the process takes to go from one such set to another are much larger 
than the time it takes to "look like" being in equilibrium with respect to the conditional distribution 
Qiv('|'S'Ar,i). Some reflection shows that this statement has considerable difficulties and cannot be 
interpreted literally, and that a precise definition of metastability is a rather tricky business (sec e.g. 
the recent discussion in [BK]). We will give a precise definition that is, however, inspired by this vague 
consideration. The main point here is that one should make precise the two time scales we alluded 
to. We will take the following attitude: to look ergodic within 3^,1, the process should have at least 
enough time to reach the "most attractive" state within Sna, while at least the times to go from 
two such states in different metastable regions should be long compared to that time. Note that this 
concept is rather flexible and allows, in general, to define metastable states corresponding to different 
time scales. 

The following definition of "stable sets" follows this ideology; however, we prefer to use certain 
probabilities rather than actual times as criteria, mainly because these are more readily computable. 
Linking them in a precise manner to times will be part of our results. We will write rf , for x G Fat, 
/ C Fjv, for the first non-zero time the process started at x arrives at I. 

Definition 2.1: A set Mn C will be called a set o/ stable points if it satisfies the following 
assumptions: There exist finite positive constants un, &jvj cjv, and tn such that they satisfy for some 
sequence en s.t. |Fjv|ejv i 0, < SNb^- 

(i) For all z e Fjv, 

P [tU. < t|] > (2.1) 

(ii) For any x ^ y G Mn, 

^[r^<r:] <a-' (2.2) 
We associate with each x G A4n its local valley 



Then 



A{x) ^ <{ z e F^ : P [r| = rX,,] = sup P [t^ = tXi„] \ (2.3) 

We will also write T^j = P[rf < t^]~^- An important characteristic of the sets / C A^jv is 
Tj = sup^g^^, Tr^.i- A simplifying assumption, that will be seen to ensure sufficient "spacing" of the 
low lying eigenvalues is that of "genericity" , defined as follows: 



^All results apply, however, also to continuous time. 
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Definition 2.2: We say that our Markov chain is generic on the level of the set A4n, if there exists 
a sequence e^r J, 0, s.t. 

(i) For all pairs x,y € Mn, o,nd any set I C AiN\{x,y} either T^j < e^Tyj or Tyj < cpfT^j. 

Each of the elements oi Ai^ in the generic case will then correspond indeed to a metastable state. 
Onr first task is to identify precisely the notion of the exit time from a metastable state. To do so, 
we define for any x & Mn the set M.n{x) = {y & Mn ■ QN{y) > Qn{x)}; these are the points that 
are even more stable than x. The metastable exit time, tx from x is then defined as the time of the 
first arrival from a; in A^jv(a;). With this notion we can formulate our main result: 
Theorem 2.3: Assume that Mn is a stable set and that the genericity assumptions are satisfied 
with Em such that rf^c^eN i 0- Set t^ = '''XiAr(x)' ^^^^ 

(i) For any x ^ Mn, 

Etx = R-^Tx^m.(.){^ + o{1)) (2.5) 
(a) For any x £ Mn, there exists an eigenvalue of 1 — Pn that satisfies 

Ax = ^(1 + 0(1)) (2.6) 

Moreover, there exists a constant c > such that for all N 

a{l - Pn)\ ^x&Mn Xx C {cbN, 1] (2.7) 

(Hi) For any x e Mn, for all t > 0, 

F[tx > tm,] = e-*(^+°(^»(l + 0(1)) (2-8) 
(iv) If ipx denotes the eigenvector of 1 — Pn corresponding to the eigenvalue Xx, then 

^.(,) = / 'I'-' < + '1^-' < ^ « (2.8) 

O(ejv), otherwise 

Remark: Explicit bounds on the error terms are given [BEGK2,EK]. 

Let us make some additional comments on this theorem. First of all, they identification of what 
constitutes a metastable exit is crucial, and, in particular the fact that these processes include the 
transition through the "saddle point", guaranteed in our case by the insistence that the process has 
actually arrived in Mn{x). Without taking this into account, the precise uncertainty principle (ii) 
could not hold. It is interesting to note that on the level of this theorem, the difficulties associated with 
the control of the passage through a saddle are not visible, and that we have the exact formula (2.5) 
for the mean exit time. Of course the difficulty is hidden in the quantities Tx^y whose computation 
is far from trivial. However, we have shown in [BEGKl] that at least in the reversible case, using 
a variational representation, very precise control of such quantities can be gained in certain setting. 
Somewhat less precise results can also be obtained in some non-reversible situations [EK] . Concerning 
our estimate on the eigenf unction, it is easy to see that [BEGK2] P[t^ < t^^^^^-^] is either very close 
to one or very close to zero, except on a set of points whose invariant measure is extremely small. 
Therefore, the corresponding right-eigenfunctions tpxi^) = Qn{z)'^x{z), are essentially proportional 
to the measure Qn conditioned on the local valley of corresponding to x (all up to errors of order ejv), 
i.e. they do indeed represent metastable measures, as suggested in [GS]. 
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3. Some ideas of the proofs. 

The first major ingredient for the proof is a representation formula for the Green's function of the 
transition matrix in terms of certain probabiUties. It implies in particular that for any / C Fjv, 

yerN\i\x ^^'^ ^ L / xj 

This formula was first derived for the reversible case in [BEGKl]. An apparently independent deriva- 
tion that also covers the non- reversible case was given rciccintly in [GM] . This formula allows to prove 
in a rather simple way, (2.5). However, the realization that this formula actually arises from a repre- 
sentation of the Greens function makes it even more useful. 

Our analysis of the spectrum of 1 — Pjv passes through the analysis of the Laplace transforms, 
Gy j{u) = Ee"*3'llT-x<T-» of transition times of process that is 'killed' upon arrival in a set / C Fjv. We 
write for the transition matrix of such a process, and we write Aj for the smallest eigenvalue of 
(1 — P^). It then turns out that all eigenvalues of (1 — Pn) below Xj can be characterised as follows: 
Set u{X) = — ln(l — A). Define the |J| x | J| matrix Gj{u) whose elements are 

5m',m - G"!j(«), m',meJ (3.2) 

Then A is an eigenvalue of (1 — Pn) below Aj, if and only if 

detGj(u(A)) = (3.3) 

This equation is rather easy to understand if | J| = 1. In this case, (3.3) becomes simply G™(m(A) = 1. 
By a simple renewal argument, one sees that G^{u) = j^^^^-j^)- Therefore, u(A) defined by (3.3) is 

the first value at which sup^g^^ G™{u) = +oo. The general formula (3.2) is somewhat less intuitive. 
Basically, one makes an ansatz for the eigenfunctions of (1 — Pn) in terms of the Laplace transforms 
of the form 

^ix) = J2 ^mG-m,A^) (3.4) 

One then finds that condition (3.3) is sufficient for the ansatz to yield eigenfunctions wit u = u(A). 
Moreover, one can show that if A is an eigenvalue, then the eigenfunctions can be represented in this 
form and (3.3) must be satisfied. 

To complete the proof one needs good control over the Laplace transforms; this is partly provided 
again by the representation of the Green's function, complemented by lower bounds on eigenvalues Aj 
obtained from a Donsker-Varadhan [DV] argument. The actual proofs are rather involved and must 
be left to the longer publications [BEGK2,EK]. Let us finally mention that the good control over the 
spectrum of (1 — Pjv) allows in turn a very good control of the analytic properties of the Laplace 
transforms which allow in turn the sharp estimates on the probability distribution of metastable 
transition times stated under (iv). 
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